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The randomly driven Burgers equation with pressure is considered as a ID model of strong tur- 
bulence of compressible fluid. It is shown that inflnitely small pressure provides a flnite effect on 
the velocity and density statistics and this case therefore is qualitatively different from turbulence 
without pressure. We establish the corresponding operator product expansion and predict the in- 
termittent velocity-difference and mass-difference PDFs. We then apply the developed methods to 
the statistics of a passive scalar advected by the Burgers field. 
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The Burgers equation with a random external force is 
considered to be the first exactly solvable model of ID 
turbulence and has been extensively studied in recent 
years Though rather simplified, this model 

can serve as a test model for some general ideas within 
the theory of strong turbulence. Recently, methods of 
quantum field theory were applied to this problem by 
A. Polyakov Q], which enabled a qualitative explanation 
of the numerical findings by A. Chekhlov and V. Yakhot 
1^. In 1^ it was shown that the Polyakov approach al- 
lows, in fact, to obtain the quantitatively correct results. 
We believe that the operator product expansion, intro- 
duced in 1^] to take into account the viscous term, is an 
adequate language to treat compressible turbulence (even 
in higher dimensions, where shock structures and associ- 
ated local dissipation persist). In the present paper we 
extend these ideas to Burgers turbulence with pressure, 
a more realistic model of ID turbulence. 

The basic equations we will study are the following: 



2 7 — 2 
XX ^ P Px 



Ut + UUx = VU 

Pt + {pu)x = , 



/ 



(1) 



where we assume pressure of the form p ^ c^p^ . It is 
the simplest model where we take into account the back 
reaction of the density field to the velocity dynamics. 
For 7 = 1, c coincides with the sound velocity. The force 
is chosen to be Gaussian with zero mean and white in 
time variance: 



{f{x,t)f{x',t'))=5{t-t')K{x-x') 
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where the k function is concentrated at some large 
scale L. We assume that steady states for both veloc- 
ity and velocity difference exist; for this we can require, 
for example, that periodic boundary conditions on a scale 
much larger that L are imposed, and that the zero har- 
monic in the n function is absent. 

The considered model has many applications for fluid 
and plasma problems; instead of kinematic pressure it 
can contain magnetic pressure and thus describe MHD 
turbulence, etc. The resulting equations for the ID case 
are quite similar to system (|l|). Their derivation and un- 
derlying physics can be found in 0,0]. These equations 



with white in time and space random forcing were used 
to analyze the MHD turbulence in ||^ . 

In this paper we appeal to the results obtained for 
Burgers turbulence without pressure in In partic- 

ular, we are interested in the velocity-difference proba- 
bility density function (PDF), where both velocities are 
taken at the same time. The physical picture, presented 
in these papers, allows us to consider such a general phe- 
nomenon as intermittency on a rigorous basis; it is related 
to the spontaneous breakdown of the Galilean invariance 
of the forced equation and to the algebraic decay of the 
PDFs. We will not repeat these arguments here, instead, 
we will concentrate on the main ideas which allow us to 
include the pressure effects into the picture. 

First, we would like to note that equation (^ devel- 
ops singularities (shocks), and that both dissipation and 
pressure terms play a regularizing role. The viscous term 
regularizes the theory at small scales and also provides 
energy dissipation. Besides this, it does not make much 
physical sense. Strong shocks have the width of the or- 
der of the mean free path, where the gradient expansion 
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used to derive the viscous term in (1^ does 



not hold. In principle, this term can be chosen in many 
different ways and the steady state can be sensitive to 
this choice. The general results we disscuss below will 
account for such a possibility. The same is true for the 
pressure term in (|l|), the role of which will be clarified 
later. 

We start by neglecting the external force and consid- 
ering the special stationary solutions of these equations: 
shock waves or kinks. Due to the Galilean invariance 
of (|l|), we can work in the system, where the shock front 
is at rest. The second equation gives p = B /u with some 
constant B. Let ui denote the fluid velocity to the left 
of the shock and U2 to the right. Obviously, there can be 
either mi > 0, U2 > or ui < 0, M2 < 0, due to continu- 
ity. Assume first that ui > 0, M2 > 0, then B > Q. For 
simplicity, choose 7 = 1. The solution for the resulting 
equation has the form: 
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where U1U2 = <? ■ From this solution one can see that 
u\ > U2 and pi < p2, so we have a kink ("— tanh"- 
Hke jump) for the velocity and an antikink ( "+ tanh"-like 
jump) for the density. For the other case, ui < 0, it2 < 0, 
we will have kinks for both velocity and density. 

Two different kinds of shocks are thus possible for the 
flux with pressure. For both of them the velocity has 
a "— tanh"-like jumps, as in turbulence without pres- 
sure, while the corresponding density profiles are dif- 
ferent: both "— tanh" and "-|-tanh" jumps are allowed. 
Physically, this corresponds to shock waves moving in 
opposite directions. This picture does not change quali- 
tatively for 7 > 1. 

Now we turn on the force and consider the dynamics 
of the shocks with small pressure, (^p^f^^i-^ ^ 1- When the 
shock is just created or two shocks with density kink and 
antikink collide, a spike appears in the density field, a 
(^-singularity, which is regularized by the pressure. One 
can say that while the viscosity regularizes the velocity 
field, the pressure term regularizes the density field. The 
presence of pressure (arbitrarily small) can then result 
in a repulsion of two colliding shocks. To describe this 
process we note that there exists a length, characteristic 
to such a collision: 



Here D is the anomaly term coming from both dissipa- 
tion and pressure in the limit of small v and c^. Without 
further assumptions, this term can not be rewritten in 
terms of the Z function. To find its structure we con- 
sider the products of operators, entering D: 
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The idea now is to consider both of them using the 
OPE. To do this, we first note that the OPE should be 
invariant under the rescaling of p: p ap. This rescal- 
ing leads simply to the rescaling of the cutoff parame- 
ter c^, and therefore should not affect the statistics at 
large scales. We have the following UV finite operators 
at our disposal: e^", Ux, p- The only Galilean invariant 
and finite combination satisfying the rescaling condition 
is pe^". Therefore, we assume the following fusion rules 
for v and c 



0: 
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If the universal scaling invariant solution of (g) exists, 
its general form is given by: 



Z = Fiy)Zoipy), F{y)^y- 
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In 
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This is a characteristic distance at which the collision 
takes place. It depends on the local velocity u, but for 
estimates one can take u ^ u-rmsi 'Urms — {LK{0)y^^ . 
For the marginal case 7 = 1, this length vanishes, and 
the collision occurs at the dissipative scale. (It is the 
case when the small pressure does not affect the statis- 
tics of velocity and can be considered as a perturbation.) 
Analogously, we introduce the dissipative length, which 
for the second order dissipative regularization takes the 
form li, ^ vju. We will assume that 7 > 1, and set v 
and (? to be small. More precisely, we will be interested 
in the field correlators at scales much larger than the 
dissipative and pressure ones. 

To proceed on the quantitative level, we introduce the 
following Z function: 

Z{p,y) - (p(a; + y/2)p(x-2//2)e^["(-+^/2)-«(--^/2)l). (5) 

As it was shown in |l[], this function satisfies the dif- 
ferential master equation, which is closed without dis- 
sipative and pressure terms. Our goal is to extend the 
assumptions on the operator product expansion (OPE), 
Galilean and scaling invariance introduced in |^ to take 
into accout these terms. Nevertheless, we would like to 
stress that all that will follow is a self-consistent conjec- 
ture. 

If the external force has a correlator niy) = 1 — 
master equation takes the form: 



the 



S^Z 
djidy 



p^y^Z = D 



(6) 



where S is some scaling exponent. Obviously, it should 
depend on the order in which we take the non-commuting 
limits I' ^ and ^ 0, since these two limits corre- 
spond to different dynamical pictures. Substituting this 
into (0) we obtain the following Eq. for Zq: 
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which formally coincides with the Polyakov equation 
for the velocity-difference PDFs of Burgers turbulence 
0. Since the two-point density correlator exists, Zo{py) 
tends to a constant as /i ^ 0. Therefore we should look 
for the solution for Zq among the family of positive, fi- 
nite, and normalizable solutions of (|To|). Such solutions 
exist for one parametric family of parameters a and S 
((5 > 0), and were obtained in 

With our present understanding we can not find the 
anomaly a from the general consideration. In principle, 
it can depend on the parameter 7 in the pressure term, 
exactly in the same manner as the form of the viscous reg- 
ularization can affect the steady state. We encounter the 
same problem in the equation for the velocity-difference 
PDF. The operator product {i^Uxx — p''~^ Px)e^^ enter- 
ing (1) should be expandable in e'^" and Uxe^^ by the 
same rescaling arguments as in (^, though the coeffi- 
cients of the expansions are undetermined. 

Nevertheless, simple arguments allow us to predict a 
rather peculiar structure of the velocity-difference PDF, 
in particular its left tail, which is determined by shocks 
This tail has an algebraic decay for Au —00, 
and is defined by the possibility of having a shock be- 
tween two points X — y/2 and x-\-y/2. To describe it we 
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first note that the dissipative and pressure lengths, as the 
functions of local velocity, are equal at some cross-over 
point M*. If y/L <C u*/urms ^ 1, there will be two differ- 
ent decay regimes: for Am 3> u*, and for Au <^ u^:. For 
example, for the quadratic viscous regularization we have 

= :/(7-l)/(7-3)/c2/(7-3)^ < 3^ the Au » 

tail will be the same as without pressure while 
the part Au ^ will be determined by the pressure. 
Moreover, this latter part will decay faster than the for- 
mer one, since the pressure term splits large shocks into 
smaller ones. This function for the regular external force 
K{y) = 1 — and for 7 < 3 is schematically presented 
on Fig. 1. 

P(u) 




M» U 

FIG. 1. Velocity-difFerence PDF for turbulence 
with pressure (7 < 3). We refer to Am simply as u. 



To complete the picture we now need to describe the 
statistics of the density field. We will describe it by in- 
troducing the probability of having mass Am between 
two points at the distance y <^ L. This PDF is defined 
for Am > 0. Presence of pressure leads to large mass 
correlation at large distances, and to the existence of the 
universal part of the mass-difference PDF. Due to the 
scaling invariance the general form of the universal part 
is expected to be P{Am,y) = p{Am/y^')/y'^, where the 
exponent r] depends on 7. Consistence with (^) fixes this 
exponent to be 77 = 1 — (5(7). We suggest that the uni- 



versal part of this PDF has an algebraic tail 

Pi^m) ^ py ^Am^pLo. (11) 
(Am)" 

where the 6 exponent is determined by the pressure reg- 
ularization. The scaling invariance breaks for large Am, 
since the total mass is conserved {Lq is the dimension of 
the system), and also for small Am, due to the existence 
of p. 

This physical picture gives a rather simple prediction 
for the mass-difference structure functions ((Am)"). Be- 
ginning with some n, these functions are non-universal. 
Their non-universal behavior corresponds to large mass 
accumulations between two points, and therefore co- 
incides with the behaviour without pressure, when all 
the mass is accumulated on shocks. The same non- 
universality should show up if we fix n and change 7. 
Numerical data for the {p{0)p{y)) correlator do reveal 
the existence of such a nonuniversality, and show that 
the critical 7 for this case is 70 « 2. For 7 < 2, the 
second order density statistics do not feel the pressure 
anomaly and F{y) = {p{0)p{y)) approaches F{y) ~ 
as we decrease c^, as it was without pressure^. Since large 
mass accumulations occur at the moments of collisions of 
shocks, one can say that the shock gas becomes strongly 
coUisional for 7 < 2. For 7 > 2, the main contribution to 
the density correlator comes from the universal part of 
the PDF ([ll|), and the corresponding exponent depends 
on the pressure regularization. Our numerical simula- 
tions seem to be in good agreement with this picture. 
We have considered the limit y ^ Ip ^ 1^, and checked 
that for 7 < 2, the exponent 25 tends to 2 as we decrease 
Ip] the same results as we observed without pressure. For 
2 < 7 < 3 the exponent 2S tends to a limit which does 
depend on 7. It is interesting that this dependence is 
linear and is rather close to 2(5 = (7 -|- l)/2. 

We emphasize once again that this general consider- 
ation does not enable us to find exponents of the tails 
of velocity-difference and mass-difference PDFs, deter- 
mined by the pressure; the exact methods should be de- 
veloped for this purpose |]8|"|l2|. Nevetheless, the pre- 



^This result has a simple qualitative explanation. One can 
show that the continuity equation with the velocity taken in 
the form of the Burgers shock, it — —Uota,n]i{xUo/v), has 

a quasistationary solution p{x,t) ~ for S> a;^ ^ 



^ -2tf/o/^^ Therefore, 

1 /"""^ 1 1 
(p(O)p(a)) ~ - / IT— 7 lTXTj^^ ~ ln(a)/a^ , 

^ Q I JL LL I I tly I CL I 

where we average over the shock position x between the two 
points —a and a. The contribution comes from small enough 
shocks, with Uo ~ v/a. We did not distinguish In corrections 
in the numerical simulations. 
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sented picture is rather suggestive and allows for numer- 
ical checks. For example, it predicts the following condi- 
tional mean for large negative Au: 



{piO)piy)\Au) ^ y'" (Au 



(12) 



where b is the anomaly for the velocity-difference PDF 
The detailed numerical simulations will be carried 
out later and presented elsewhere. 

In conclusion we would like to briefly discuss another 
important application, a passive scalar advected by the 
Burgers velocity field. It can be density, concentration or 
charge of admixture particles, temperature field or weak 
magnetic field in a turbulent flow, etc. Using the meth- 
ods for obtaining the two point density correlator one 
can immediately predict the two point correlator for the 
passive scalar. Indeed, if we neglect the driving and dissi- 
pative terms, the derivative of the passive scalar field 
is governed by the continuity equation: 







(13) 



Admixture density and charge correspond to the field , 
admixture concentration, temperature field and magnetic 
field correspond to the field T . Though the problems are 
not exactly identical (the T field does not affect the ve- 
locity dynamics and should not, in general, be strictly 
positive), the scale invariant solution and the anomalies 
should be written in the same form. The basic equation 
for the passive scalar advection is: 



Tt -f uT^ = viT^j, -f /i , 



(14) 



where /i is white in time Gaussian with zero mean and 
the correlator Ki{y) = 1 — y. This correlator is assumed 
to be concentrated at some large Lt- By analogy with 
the density, one can write the closed equation for the 
following Z function: Z2(y,/i) — {Tx{x + y/2)Tx{x — 
y/2)e^'^"}. Considering the anomaly: 

we note that the translation symmetry of equation (|l^ ) 
(T — > T -f a) is spontaneously broken by the external 
force (Trms = ('^i(O)iT)^^'^), but restores for large Lt- 
We also assume the scaling symmetry of the fluctuations. 
The presence of such symmetries means that the OPE 
should be proportional to T^. Therefore, the only possi- 
ble flnite and Galilean invariant form is: TxC^'^. Then, 
we have the following equation for the Z2 function: 



dfidy 



- p^y^Z2 = aZ2 - At" (0)^2 
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where Z2 is the Z function for the velocity difference, and 
is known. Now we assume that for y there exists the 
scaling solution, Z2 = y~^^ Zq. This can only be true 
if for small y we can neglect the last term in Eq. (15): 
aZ2 ^ ^1(0)^2, i.e. the solution should be dominated by 



the zero mode of (p^). It is important that this inequal- 
ity hold uniformly for any /i. The following consideration 
shows that this can be satisfied only in the case Zq = Z2. 
Consider the asymptotics of the velocity-difference PDF 
(which is the Laplace transform of Z2): 
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One can obtain the asymptotics of the Laplace trans- 
form of Zq by simply changing b 6 in these expressions. 
Comparison of these asymptotics shows that the inequal- 
ity holds uniformly with respect to u only when S — b. 
In other words, the only zero mode surviving in the limit 
y/i <C 1 is the mode with S = b. 

One can therefore predict for the second order struc- 
ture function: 



mO)-T{y)r)^y 
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For each a in the short-distance expansion K{y) = 1 — y" 
some particular solution b = b{a) realizes. For the case 
of a regular external force, a = 2, it is argued in that 
6=1/2, and we obtain ((r(0) - T{y)f) - ?/ln(?/). 
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